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ABSTRACT 

This paper is the third in the series of papers on M81 globular clusters. In this paper, we present 
spatial and metal abundance properties of 95 M81 globular clusters, which comprise nearly half of 
all the M81 globular cluster system. These globular clusters are divided into two M81 metallicity 
groups by a KMM test. Our results show that, the metal-rich clusters did not demonstrate a centrally 
concentrated spatial distribution as ones in M31, and metal-poor clusters tend to be less spatially 
concentrated. In other words, the distribution of the metal-rich clusters in M81 is not very similar to 
that of M31. Most of the metal-rich clusters distribute at projected radii of 4-8 kpc. It is also noted 
that the metal-rich clusters distribute within the inner 20 kpc, and the metal-poor ones do out to radii 
of ~ 40 kpc. Like our Galaxy and M31, the metallicity distribution of globular clusters in M81 along 
galactocentric radius suggests that some dissipation occurred during the formation of the globular 
cluster system, i.e. smooth, pressure-supported collapse models of galaxies are unlikely to produce 
such radial distribution of metallicity presented in this paper. There is not evident correlation between 
globular cluster luminosity and metallicity in M81 globular clusters. The overwhelming conclusion of 
this paper seems to be that a more complete and thorough cluster search is needed in M81. 

Subject headings: galaxies: individual (M81) - galaxies: star clusters - globular clusters: general 



1. INTRODUCTION 

An understanding of galaxy formation and evolution is 
one of the principal goals of modern astrophysics. Glob- 
ular clusters are fossils of the earliest stages of galaxy 
formation and evolution. They are bright, easily rec- 
ognized packages containing a stellar population with a 
homogeneous abundance and age. So, their integrated 
properties of location, abundance, and kinematics pro- 
vide valua ble clues to the natu re and duration of galaxy 
formation (Bar mbv et aLll^OOOT ). 

The metallicity distribution of globular clusters is 
of particular importance in deepening our knowledge 
of the dynamical and chemical evolution of the par- 
ent galaxies. For example, the globular clusters of 
many elliptical galaxies show multi-modal metallicity 
distributions, suggesting that multiple star formation 
episo des occurred in these elliptical galaxies in the 
past (|Zepf fc Ashmanlll993l:rBarmbv et al.ll2000f) . Great 
progress has been made in the past decade in our un- 
derstanding of globular cluster systems of galaxies, es- 
pecially the discovery that many galaxies possess two or 
more distinct sub populations of globular clusters (e.g., 
IWest et"aLl 12004 and references therein). Based on 
the data from Hubble Space Tel es cope (HST) a rchive , 
iGebhardt fc Kissler-Pat^ (119991) . lLarsen et all (|200lD 
and iKundu fc Whitmord ((2001) presented that many 
large galaxies possess two or more subpopulations of 
globular clusters that have quite d i fferen t chemical com- 
positions. Recently, iPeng et all (|2006l ) presented the 
color distributions of globular cluster systems for 100 
early-type galaxies observed in the Virgo Cluster Survey 
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with the Advanced Camera for Surveys (ACS) on the 
HST, and found that, on average, galaxies at all lumi- 
nosities in their study appear to have bimodal or asym- 
metric GC color/metallicity distributions. The presence 
of color-bimodality indicates that there have been at least 
two major star-forming mechanisms in the histories of 
g alaxie s. 

ICotel (|1999f ) presented a metallicity distribution of 133 
Galactic globular clusters that apparently shows two 
peaks (i.e., two distinct metal-poor and metal-rich glob- 
ular cluster populations). A double-Gaussian can best 
fit these two subpopulations, the mean metallicity values 
are —1.59 and —0.55 dex, res pectively. Using the d ata for 
247 globular clusters in M31, lBarmbv et all (|2000D stud- 
ied the metallicity distribution, which is asymmetric, im- 
plying the possibility of bimodality. Then they applied 
a KMM algorithm show ing that the metalli city distri- 
bution is really b imodal. iPerrett et al.l (120021) confirmed 
the conclusions of lBarmbv et all (|2000f) . ]Ma et all (|2005l ) 
showed that the intrinsic B and V colors and mctallicities 
of 94 M81 globular clusters are bimodal, with metallicity 
peaks at [Fe/H] w —1.45 and —0.53, similar to what we 
find for the Milky Way and M31 globular clusters. 

M81 is one of the nearest Sa/Sb-type spiral outside the 
Local Group, very similar to M31, and roughly as mas- 
sive as the Milky Way. So, beyond the Local Group, 
it is a good candidate for reaching a detailed study 
of spiral galaxy globular cluster syste m for comparison 
to the Milky Way and M31 system. iBrodie fc Huchral 
(|1991[ ) derived spectroscopic metallicities for eight glob- 
ular clusters in M81 a nd presented the sampl e mean of 
[Fe/H] = -1.46±0.31. iPerelmuter et all (| 19951 ) obtained 
low signal-to-noise spectra of 82 candidates, 25 of which 
were confirmed as bona fide M81 globular clusters. They 
derived the mean metallicity to be [Fe/H] = — 1.48±0.19 
both from the weighted mean of the individual mctallic- 
ities, and directly from the composite spectrum of the 
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25 confirmed globular clusters. To maximize the success 
rate of the globular cluster candidate list for the ongoing 
spectroscopic observations, iPerelmuter fe Racind (1995) 
used an extensive database that included photometric, 
astrometric, and morphological information on 3774 ob- 
jects covering over a > 50 arcmin diameter field centered 
o n M81 to reveal 70 glo bular cluster candidates. 

iSchroder et al.1 (|2002f ) presented moderate-resolution 
spectr oscopy for 16 globular cluste r candidates from the 
list in IPerelmuter fc Ra cine i|1995[) . and confirmed these 
16 candidates as bona fide globular clusters. They also 
obtained metallici ties for 15 of th e 16 g lobular clusters. 
From their results, ISchroder et al.l (|2002f ) concluded that 
the M81 globular cluster system is very similar to the 
Milky Way and M31 systems, both chemically and kine- 
matically. 

With the superior resolution of the HST, M81 is close 
enough for its clusters to be easily resolved on the ba- 
sis of image structure (jChandar et alj 120011 ). Thus, us- 
ing the B, V, and I ban ds of HST Wide Fiel d Plane- 
tary Camera 2 (WFPC2). IChandar et all (|2001h imaged 
eight fields covering a total area of ~ 40 arcmin 2 , and de- 
tected 114 compact star clusters in M81, 59 of which are 
glo bular clusters. Based on the estimated intrinsic col- 
ors, IChandar et"afl (|2004h found that the M81 globular 
cluster system has an extended metallicity distribution, 
which argues the prese nce of both meta l-rich and metal- 
poor globular clusters. iMa et all (|2005f ) then confirmed 
this conclusion. 

The outline of the paper is as follows. In § 2 we provide 
some statistical relationships. The summary is presented 
in § 3. 

2. PROPERTIES OF GLOBULAR CLUSTERS IN M81 
2.1. Sample of globular clusters 

In the first paper of our series, IMa et all (|2005l) 
studied the distributions of intrinsic B and V colors 
and metallicities of 95 M81 g lobular clusters which are 
from IPerelmuter et all (|1995h . IChandar et alj (|200lD and 
ISchroder et alj (|2002f ). This cluster sample includes 
nearly half of the M81 total globular cluster popula- 
tion. About the M81 total globular cluster population, 
IPerelmuter fe Racind (|1995l ) estimated it to be 210 ± 30 
by the BVR photometric, a strometric, and i mage struc- 
ture study of the M81 field; IChandar et~aTl (|200lD indi- 
cated the total number of globular clusters in M81 to be 
211 ±29 using globular cluster estimates in various annu- 
lar bins and correcting for incompleteness. It is difficult 
to detect and confirm globular clusters beyond the Local 
Group before HST appears. In fact, even in the Local 
Group, it is also not easy to detect and confirm glob- 
ular clusters. For example, the total number of glob- 
ular clusters in M 3 1 was estimated to be 460 ± 70 by 
iBarmbv fc Huchral (|2001h . the largest number of globu- 
lar clusters used to study the metal abundance properties 
of the M 31 globular c l usters includes 301 clusters col- 
lected bv lPerrett et al.l (|2002f ). a little more than half of 
the total number. Beyond the Local Group, the glob ular 
cluster sample in M81 collected by IMa et all (|2005t ) in- 
cludes the largest number of globular clusters comparing 
to the total globular clusters in the ho st galaxy. 

In the second paper of our series, IMa et all (l2006h 
presented the spectral energy distributions of 42 M81 



Fig. 1.— The image of M81 in filter BATC07 (5785A) and the 
positions of the sample star clusters. The center of the image is 
located at RA = 01 h 33 m 50?58 Dec=30°39'08'.'4 (J2000.0). North 
is up and east is to the left. 

globular clusters selected from IMa et ail ((2005) in 13 
intermediate-band filters from 4000 to 10000A, using the 
CCD images of M81 observed as part of the Beijing- 
Arizona- Taiwan-Connecticut (hereafter BATC) multi- 
co lor survey of the sky, and confirmed the conclusions 
of ISchroder et all |2002) that, M81 contains clusters as 
young as a few Gyrs, which were also observed in both 
M31 and M33. In this paper, we will study the spa- 
tial and metal abundance properties of the M8 1 globular 
cluste rs using the sample globular clusters of IMa et"ai1 
(2005!) ■ Figure 1 is the image of M81 in filter BATC07 
(5785A) of BATC multicolor survey of the sky, the circles 
indicate the positions o f the sample clu sters. 

As mentioned above. IMa et al.l (|2005l ) studied the dis- 
tributions of intrinsic B and V colors and metallicities 
of these 95 M81 globular clusters, and first found that 
the abundance distribution of the globular cluster sys- 
tem is consistent with a bimodal distribution with peaks 
at [Fe/H] -1.45 and -0 53 b ased on the KMM al- 
gorithm of I Ashman et al.l (|1994l ). It is true that, the 
appearance of the histogram can be ambiguous and mis- 
le ading with binned dat a, however, the KMM algorithm 
of lAshman et alJ ([T994) is a robust meth od of analysis 
witho ut relying on binning methods (see iPerrett et al.1 
2002). KMM mixture modelling works under the as- 
sumption that the sample data are independently drawn 
from a parent populatio n that comprise s a mixture of N 
Gaussian distributions. IMa et all (|2005l ) presented that, 
for M81 globular clusters, the posteriori probabilities of 
group membership returned by the KMM algorithm as- 
signed 74 clusters to the metal-poor population and 20 
to the metal- rich population distribution 4 . 

2.2. Spatial distribution 

Figure 2 shows the projected spatial distributions 
of the metal-poor and metal-rich globular clusters in 
M81. The distance modulus for M81 is adopted to be 
27.8 (|Freedman et al.l Il99l IChandar et~aH l200ll) . We 
adopted the incl ination and pos i tion a ngles to be 59° and 
157° of M81 as IChandar et all (|200lD did, respectively. 
When the line of intersection (i.e. the major axis of the 
image) between the galactic plane and tangent plane is 
taken as the polar axis, it is easily proved that: 



r = pJl + tan 2 7 sin 2 9 (1) 

and 



where r and <p are the polar co-ordinates in the galactic 
plane, and p and 9 are the corresponding co-ordinates 
in the tangent plane, and 7 is the inclination angle of 
the galactic disk. Using formula (1), we can obtain the 
distances of our sample clusters from the center of M81, 

4 Since the globular cluster 96 of Chandar ct al. (200l|) has very 
high (B — V)o ((B — V)o = 1.778), and the metallicity ob tained us- 
ing th e color-metallicity correlation is too rich (0.95 dex). lMa et al.l 
(2005) do not include it when performing the KMM test. 
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which are listed in Table 1. In Table 1, we also listed th e 
metallicities of the sample clusters from lMa et al.l (|2005f ). 
From Figure 2, it is clear that the metal-rich globular 
clusters in M81 are not as centrally conce ntrated as the 
meta l -rich globular clust ers of M31 are (|Huchra et al.l 
11991b IPerrett et aD l2002f ). Figure 3 presents the his- 
togram for the metal-poor and metal-rich globular clus- 
ters in M81. It shows that most of the metal-rich clusters 
distribute at projected radii of 4-8 kpc. It is also noted 
that the metal-rich clusters distribute within the inner 20 
kpc, and the metal-poor ones do out to radii of ~ 40 kpc. 
In the Milky Way, the metal-rich GCs reveal significant 
rotations and have historically been associated w ith the 
thick-disk system (IZinnl fl985t lArmandrofil Il989f); how- 
ever othe r works flFrenk fc White! Il982t iMinnitil Il995t 
lC6telll999t iForbes et al.ll200lD suggested that metal-rich 
GCs within ~ 5 kpc of the Milky Way Galactic cen- 
ter a re better associated with the bulge and bar. In 
M31. lElson fc Walterboi (|1988f ) showed that the metal- 
rich clusters constitute a more highly flatted system than 
the me tal-poor ones, and a ppear to have disklike kine- 
matics; iHuchra et alj (|1991f ) showed that the metal-rich 
GCs are pr eferentially close to the galaxy center. At the 
same time. lHuchra et akl (|1991l) showed that the distinc- 
tion between the rotation of the metal-rich and metal- 
poor clusters is mo st apparent in the inner 2 kpc. So, 
IHuchra et al.l (|199lD concluded that the rich-metal clus- 
ters in M31 appear to form a central rotatin g disk system. 
With the largest sample of 321 velocities. IPerrett et al.l 
(2002) provided a more comprehensive in yestigation on 
the ki nematics of the M31 cluster system. IPerrett et al.l 
(2002) showed that, the metal-rich globular clusters of 
M31 appear to constitute a distinct kinematic subsys- 
tem that demonstrates a centrally spatial distribution 
with a high rotation amplitude, but does not appear sig- 
nificantly flattened, which is consistent wi th a bulge pop- 
ulatio n. It is of interests to mention that, Schroder et alJ 
(2002) performed a maximum-likelihood kinematic anal- 
ysis on 166 M31 clusters of iBarmbv et al.1 (|2000D and 
found that the most significant difference between the 
rotation of the metal-rich and metal-poor clusters occurs 
at int ermediate proj e cted galactocentric radii. Espe- 
cially, ISchroder et al.l (|2002l ) presented a potential thick- 
disk population among M31's m etal-rich globular clus - 
ters. For M81 globular clusters, ISchroder etTaTl (|2002ft 
performed a kinematic analysis of the velocities of 44 
M81 globular clusters, and strongly suggested that the 
metal-rich clusters are rotating in the s ame s ense as the 
gas in the disk of M81. ISchroder etTaTl (|2002f ) concluded 
that, although their cluster sample is not large enough to 
make a direct comparison between metal-rich and metal- 
poor clusters in specific radius ranges, the conclusion 
with M81's metal-rich clusters at intermediate projected 
radii being associated with a thick disk in M81 is cor- 
rect. It is true that, from Figure 3, most of metal-rich 
clusters distribute at projected radii of 4-8 kpc. So, at 
least, we can conclude that most of the metal-rich clus- 
ters in our sample are not associated with a bulge cluster 
system of M81; they may associate with a thick disk in 
M81 as indicated by ISchroder et all (|2002f h The sample 
clusters of this paper include all the sample clusters of 
ISchroder et al.l (12002!) . But, except for the sample clus- 
ters of ISchroder et al.l (|2002f ) , the other clusters have no 
published radial velocity estimates. So, obviously, more 
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Fig. 2. — Spatial distributions of the metal-rich and metal-poor 
globular clusters. 
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Fig. 3. — Radial distributions of the metal-rich and metal-poor 
globular clusters. 

kinematic and metallicity data are needed for globular 
clusters in M81 to determine if the inner metal-rich GCs 
have kinematic properties that are consistent with the 
bulge and the metal-rich GCs at projected radii of 4-8 
kpc are associated with a thick disk in M81. 

2.3. Metallicity gradient 

The presence or absence of a radial trend in globu- 
lar cluster metallic ities is an importan t test of galaxy 
formation theory (Barmbv et al. 2000). If a galaxy 
forms as a consequence of a monolithic dissipative and 
rapid collapse of a single massive, nearly-spherical spin- 
ning gas cloud in which the enrichment timescale is 
shorter than the collapse time, the halo stars and glob- 
ular clusters should show large-scale metall icity gradi- 
ents (lEggen et all [ 1962: B armbv et al.l 12000): however, 
iSearle fc Zinnl (|1978h presented a chaotic scheme in the 
early evolution of a galaxy, in which loosely bound pre- 
enriched fragments merge with the main body of the 
proto-galaxy over a significant period, so there should 
be ho mogeneous metallic ity distribution. For the Milky 
Way, lArmandrofil (|1989t ) showed some evidence that 
metallicity gradients with both distance from the Galac- 
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Fig. 4. — Metallicity as function of projected radius for M81 glob- 
ular clusters. Black circles indicate the clusters with spectroscopic 
mctallicities with uncertainties smaller than 1.0 dex. 

tic plane and distance from the Galactic center were 
present in the disk cluster system. For M31, there are 
some inconsistent conclusions, such as Ivan den Bergh I 
(1969) showed that there is little or no evidence for 
a correlation between metallicity and projected radius, 
but most of his clu sters were inside 50'; however , 
some authors (see e.g. iHuchra et alJll982t ISnajovlll988l 
iHuchra et all 119911 : iPerrett et alj l2002f) presented that 
there is evidence for a weak but measura ble metallicity 
gradie nt as a function of projected radius. iBarmby et al.l 
(2000) confirmed the latter result based on their large 
sample of spectral metallicity and color-derived metal- 
licity. Figure 4 plots the metallicity of the M81 glob- 
ular clusters as a function of galactocentric radius, in 
which black circles indicate the clusters with spectro- 
scopic metallicities with uncertainties smaller than 1.0 
dex. Clearly, the dominant feature of this diagram is the 
scatter in metallity at any radius. At the same time, our 
sample clusters are mainly distributed in the inner 10 
kpc. So, it is difficult to determine the metallicity gra- 
dient. It is true that, smooth, pressure-supported col- 
lapse models of galaxies are unlikely to produce a result 
like this. However, in order to present a quantitative 
conclusion, we made least-squares fits: the total sample 
of globular clusters does not have a significant metallic- 
ity gradient (-0.009 ± 0.009 dex kpc -1 ), and the clus- 
ters with spectroscopic metallicities with uncertainties 
smaller than 1.0 dex have a marginally significant gradi- 
ent (—0.01 8 ± 0.01 dex kpc" 1 ). This result is in agree- 
ment with iKong et al.l (|2000f ). who obtained metallicity 
maps of M81 field by com paring simple stellar popu lation 
synthesis models of BC96 ()Bruzual fc Cha riot 1996) with 
the integrated photometric measurements of the BATC 
photometric system, and did not find, within their errors, 
any obvious metallicity gradient from the central region 
to the bulge and disk of M81. But, we should emphasize 
that, in the least-squares fits of this paper, the metal-rich 
clusters seem to act an important part in determining 
the metallicity gradient. From Figure 4, we can also see 
the decre ase in the "upper en velo pe" of metallicity re - 
ported bv IHuchra etaU (j 199 If ) and lBarmbv et all (|2000l ) 
for M31 globular clusters. In fact, we do a least-squares 
fits for the metal-poor clusters of M81, the metallicity 
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Fig. 5. — Metallicity versus absolute magnitude for M81 globu- 
lar clusters. Black circles indicate the clusters with spectroscopic 
metallicities with uncertainties smaller than 1.0 dex. 

gradient is —0.006 ± 0.006 dex kpc -1 . This metallic- 
ity gradient may not have any statistical meanings. It 
is clear that, the small sample of rich- metal clusters in 
M81 cannot present any firm conclusion about metallic- 
ity gradient. 

2.4. Metallicity versus absolute magnitude 

The correlation between cluster mass (or luminosity) 
and metallicity is important in globular cluster formation 
theory. It is generally believed that, if self-enrichment is 
important in globular clusters, the most massive clusters 
could retain their metal-enriched supernova ejecta, so the 
metal abundance should increase with cluster mass; the 
opposite is true if cooling from metals determines the 
temp erature in the cluster-forming clouds ( jBarmbv et al.l 
2000). The self-enrichment of GCs has been in detail 
studi ed in some aspects (see detail from IStrader et al.1 
2006). However, it is interesting of mentioning the 
model of globula r clust er self-enrichment developed by 
iParmentier et all (jj999). In this model, cold and dense 
clouds embedded in the hot protogalactic medium are 
assumed to be the progenitors o f galactic halo globular 
cluste rs. Based on this model, IParmentier &: Gilmord 
(|2001h presented that the most metal-rich proto-globular 
clusters are the least massive ones. 

For M31 globular clusters, IHuchra et ail (|1991f) first 
presented the metallicity versus apparent magnitude for 
150 M31 clusters, and did not find any trend of metallic- 
ity with luminosity; then, IBarmby et alj (|2000f ) showed 
metallicity versus dereddened apparent magnitude using 
a large cluster sample i ncluding 247 object s, and con- 
firmed the conclusion of IHuchra et al.1 (|1991f ). 

Figure 5 shows metallicity versus absolute magnitude 
for M81 globular clusters, in which black circles indi- 
cate the clusters with spectroscopic metallicities with un- 
certainties smaller than 1.0 dex. The metallicity and 
absolute magni tudes for the sample clusters are from 
iMa et all (|2005f ) . It is true that there is not obvious trend 
of metallicity with luminosity as M31 GCs do. Least- 
squares fits show no evidence for a relationship between 
luminosity and metallicity in M81 clusters. 

As we know, HST provides a unique tool for study- 
ing globular clusters in extragalaxies. Recently, based 
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on the ACS on t he HST, lHarris et all (|2006t ) and 
iStrader et al.1 (|2006[ ) found that, in giant ellipticals such 
as M87, NGC4649 and NGC7094, luminous blue GCs 
(i.e. metal-poor GCs) reveal a trend of having redder 
colors, such that more massive GCs are more red (metal- 
ric h). This trend is refer red to as a 'blue tilt' (also 
see iBrodie fc Strader|[20M) . This 'blue tilt' was inter - 
preted as a result of self-enrichment (jStrader et al.ll2006h . 
IStrader et al.1 (|2006D speculatively suggest ed that these 
GCs o nce possessed dark matter halos. ISpitler et al.l 
(2006) subsequently found that this 'blue tilt' is also 
true in the Sombrero spiral galaxy (NGC4594) and may 
extend to less luminous GCs with a somewhat shal- 
lower slope that was deriv ed by lHarris et alj (l2006f) and 
IStrader et all (|2006f) . As ISpitler et all (|2006l ) pointed 
out that, the Sombrero provides the first example of 
this trend in a spiral galaxy and in a galaxy found in 
a low-density galaxy environment. However, in these 
ACS studies, the metal-rich (r ed) GCs did not s how this 
corresponding trend (also see lBekki et al.ll2007f ). Based 
on high-reso lution cosmo l ogical simulation with globu- 
lar clusters, iBekki et al.1 (|2007f ) investigated formation 
processes and physical properties of globular cluster sys- 
tem in galaxies, and found that, luminous metal-poor 
clusters would emerge a correlation between luminosity 
and metallicity if they originated from nuclei of low- 
mass galaxies at high z. In fact, in the simulations of 
IBekki et al.l (|2007f ). the 'simulated blue tilts' emerge from 
the assumption that luminous metal-poor clusters origi- 
nate from stellar galactic nuclei of the more massive nu- 
cleated galaxies with a lu minosity-me t allicity relation. 
So, it is evident that, in Bek ki et ail {2007;), galaxies 
which experienced more accretion/merging events of nu- 
cleated low-mass galaxies are more likely to show a blue 



tilt (see details from IBekki et al.ll2007f) . 

3. SUMMARY 

In this paper we present spatial and metal abundance 
propertie s of 95 M81 glo bular clusters, which are col- 
lected by iMa et al.l (|2005l ). This cluster sample includes 
nearly half of the M81 total globular cluster population, 
is the largest one comparing to the total globular clus- 
ters of the host galaxy beyond the Local Group. Our 
conclusions are as follows: 

1. The metal-rich clusters did not demonstrate a cen- 
trally concentrated spatial distribution as ones in M31, 
and metal-poor clusters tend to be less spatially con- 
centrated. Most of metal-rich clusters distribute at pro- 
jected radii of 4-8 kpc. We can conclude that most of 
the metal-rich clusters in our sample are not associated 
with a bulge cluster system of M81; the y may associate 
with a thick disk in M81 as indicated by Schrode r et al.l 
(l27)02h . 

2. The globular clusters in M81 have a small radial 
metallicity gradient like M31 and our Galaxy, suggesting 
that some dissipation occurred during the formation of 
the globular cluster system. 

3. There is not obvious trend of metallicity with lumi- 
nosity in M81 globular clusters. 
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tional Basic Research Program of China (973 Program) 
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